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ABSTRACT 

^ The interaction of a relativistic fireball with its ambient medium is described 

> ; 

, through two shocks: a reverse shock that propagates into the fireball, and a 

forward shock that propagates into the medium. The observed optical flash of 

(•~^ ' GRB 990123 has been considered to be the emission from such a reverse shock. 

I The observational properties of afterglows suggest that the progenitors of some 

, GRBs may be massive stars and their surrounding media may be stellar winds. 

^ '■ We here study very early afterglows from the reverse and forward shocks in 

O 

^ ' winds. An optical flash mainly arises from the relativistic reverse shock while 



a radio flare is produced by the forward shock. The peak flux densities of 
optical flashes are larger than 1 Jy for typical parameters, if we do not take 

•i-H . 

I into account some appropriate dust obscuration along the line of sight. The 

5h ' 

' radio flare always has a long lasting constant flux, which will not be covered 

up by interstellar scintillation. The non-detections of optical flashes brighter 
than about 9th magnitude may constrain the GRBs isotropic energies to be 
no more than a few lO'^^ ergs and wind intensities to be relatively weak. 

Key words: gamma-rays: bursts — hydrodynamics — relativity — shock 
waves 



1 INTRODUCTION 

Gamma-ray bursts (GRBs) have been profoundly understood since the discovery of after- 
glows in 1997 (Wijers 1997; Piran 1999; van Paradijs, Kouveliotou & Wijers 2000; Meszaros 
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2002). The hydrodynamic evolution of a GRB remnant is well described by an ultra- 
relativistic forward shock that sweeps into an interstellar medium (ISM) and slows down to 
the non-relativistic phase while the resulting afterglow is due to synchrotron radiation of 
electrons accelerated by the shock. 

However, the emission from a reverse shock propagating into the shell of the GRB ejecta, 
was also predicted (Meszaros & Rees 1997; Sari & Piran 1999a). Before the reverse shock 
crosses the shell, the shocked shell matter carries an amount of energy comparable to that of 
the forward shocked interstellar medium. The prompt optical flash of GRB 990123 (Akerlof 
et al. 1999) motivates investigations of emission from the reverse shock (Meszaros & Rees 
1999; Sari & Piran 1999b; Kobayashi & Sari 2000; Kobayashi 2000). The hght curve of an 
optical flash in a uniform ISM has been well studied in both thin shell and thick shell cases 
(Kobayashi 2000) and it is found that the optical flash of GRB 990123 is due to the radiation 
from a mildly relativistic reverse shock (Kobayashi & Sari 2000). Recently, the contribution 
of reverse shocks to early radio afterglows was also discussed by Soderberg & Ramirez-Ruiz 
(2002). 

There is increasing evidence for the association between some GRBs and type Ib/c 
supernovae. The most direct evidence is the association of GRB 980425 with SN 1998bw 
(Galama et al. 1998). This association is strengthed by its radio emission after the burst (Li 
& Chevalier 1999). Up to date several GRB afterglows, e.g., GRBs 970228 (Reichart 1999), 
970508 (Sokolov 2001), 980326 (Bloom ct al 1999), 990712 (Sahu et al. 2000), 991208 (Gastro- 
Tirado et al. 2001), 000911 (Lazzati et al. 2001), and 011211 (Bloom et al. 2002), show a 
supernova component in their light curves peaking ai — {1 + 2;)tsN,peak, with tsN,peak ~ 15 
days for SN Ic. One of the most popular central engines of GRBs now is beheved to be the 
collapse of massive stars (Woosley 1993; Paczyhski 1998). This coUapsar/hypernova scenario 
naturally explains the broad and shifted iron-group emission lines in GRBs (McLaughlin, 
Wijcrs & Brown 2002). The iron line emission can be attributed in the envelope models 
either from the irradiated shallower layers on the funnel wall along the rotation axis by a 
continuing but decaying relativistic outflow (Rees & Meszaros 2000) or from the dissipation 
of the magnetized thermal bubble expanding into H envelope of the progenitor, which results 
from the interaction between the jet with continuing energy injection and the stellar envelope 
(Meszaros & Rees 2001). The place of the iron line emission in the envelope models is at H 
envelope radius r ~ 10^^ cm, where a small amount of iron, ~ 10~® Mq and ~ 10~^ M©, can 
meet the requirement. Other iron line emission models require a large line emission radius, 
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r ~ 10^^ cm, and much larger iron mass > 0.06 Mq (Vietri et al. 2001; Bottcher & Fryer 
2001). The measured emission lines in the X-ray afterglow of GRB 011211 favors the latter 
(Reeves et al. 2002). The pre-explosion wind structure can survive in envelope models but 
may be destroyed in the other models due to either the ejection of SN prior to GRB or the 
common envelope evolution in the binary system. Detailed calculations and more sensitive 
data on X-ray afterglow spectral features will distinguish between the various line emission 
models (Ballantyne & Ramirez- Ruiz 2001). 

The light curve steepening of some well-observed GRB afterglows deviates from the 
predictions of the standard model (Sari, Piran & Narayan 1998). This deviation leads to 
discussions of the environment and geometry effects. Previous studies have discussed the 
afterglow from a spherical blast wave in a wind (Dai & Lu 1998; Chevaher & Li 1999, 
2000) and from a jet in a wind (Livio & Waxman 2000; Dai & Gou 2001; Gou et al. 2001). 
The wind environments of GRBs 970228, 970508, 991216 and 000301c arc recognized in the 
spherical-wind model (Chevalier & Li 2000; Li & Chevalier 2001). In the jet model, ten 
broadband GRB afterglows (i.e., GRBs 970508, 980519, 990123, 990510, 991208, 991216, 
000301c, 000418, 000926, and 010222) have been well fitted (Panaitescu 2001; Panaitescu & 
Kumar 2001, 2002). Among them, only GRBs 980519, 990123 and 990510 are inconsistent 
with wind environments. However, we also noticed that Frail et al. (2000) gave a better fit for 
their first 63 days observations of the radio afterglow of GRB 980519 in the spherical-wind 
model. 

Information on the surrounding environment of a GRB can also be obtained by consider- 
ing the evolution of the soft X-ray absorption property according to the surrounding density 
distribution (Lazzati & Perna 2002). The expected absorption of a wind environment will 
quickly drop to the undetectable level in no more than 1 second. It conflicts with the evi- 
dence of variable soft X-ray absorption in three GRBs detected by 5e^^oSAX. Nonetheless, 
the wind environments can not be excluded by this method due to the limited samples. 
It may give some evidence for another type of GRBs environments, such as dense clouds. 
Superbubbles which are proposed as the sites of GRBs have been expected to explain the 
obtained low densities n ~ 10~^ cm~^ (Scalo & Wheeler 2001). It requires a very weak 
stellar wind at least 6 orders weaker than that of the progenitor of SN 1998bw to insure 
non-contamination of the supperbubble by the wind at a typical afterglow radius. The wind 
profile n oc may be terminated at some radius r ~ 10^^ cm, where the wind sweeps up a 
mass of ISM comparable to that of the wind and leads to a quasi-homogeneous dense shell 
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with the mean density ~ 10^ cm~^ - 10^ cm~^ and the thickness ~ 10^'' cm (Ramirez-Ruiz 
et aL 2001). This termination of the wind profile will not affect our analysis of the dynamic 
evolution of very early afterglows, but has the potential reduction of the optical flash. In 
this paper we perform a more detailed calculation of the very early afterglows of GRBs in 
wind environments including the reverse shock effect. In §2 we give the hydrodynamics of 
forward and reverse shocks. In §3 early afterglow light curves at optical and radio bands are 
calculated. We discuss some observational implications of our results in §4. The synchrotron 
self- absorption of electrons with a broken power law energy distribution is discussed in the 
appendix. 

2 HYDRODYNAMICS OF FORWARD AND REVERSE SHOCKS 

We consider a uniform and cold relativistic coasting shell with rest mass Mq, the energy 
E, the initial Lorentz factor rj — E/Mqc^, and the observed width A, which sweeps up 
a wind with the particle number density rii — Ar~^ {k — 2), where A — M/AirmpVyj — 

3 X lO^M^cm-i and = (M/1 x 10-^MQyr-^){v^/10\m/s)-^ (Chevalier & Li 1999, 
2000). Two shocks develop: a forward shock propagating into the wind and a reverse shock 
propagating into the shell. There arc four regions separated by the two shocks: (1) the 
unshocked wind, (2) the shocked wind, (3) the shocked shell material, (4) the unshocked 
shell material. Prom the shock jump conditions and equality of pressures and velocities along 
the contact discontinuity, the Lorentz factor 7, the pressure p, and the number density n 
in both shocked regions can be determined by rii, 714, and rj (Blandford & McKee 1976, 
hereafter BM). 

We define the Sedov radius, at which the swept-up wind's rest-mass energy equals the 
initial energy E of the shell, i.e., / — MqT]/ {i.T:Amp). Assuming that the shell doesn't enter 
the spreading stage, the width A is a constant, and ^4 = Mo/(47rmpr^A77) oc r~^. The 
property of the reverse shock is largely determined by a parameter defined as / = n^/rii 
(Sari & Piran 1995), which can be further given by 

As shown by Sari & Piran (1995), if / » r/^, or the baryon loading Mq :$> AnnipAArj^, 

the reverse shock will be a Newtonian reverse shock (NRS), and the reverse shock will be a 
relativistic reverse shock (RRS) ii f ^rf or Mq <^ -iT^mpAAyf . For RRS, the Lorentz factor 
of the shocked shell relative to the unshocked shell is 73 = 7y^/^/~^/^/-\/2, and the Lorentz 

© 0000 RAS, MNRAS 000, 000-000 



factors of shocked matter are 72 = 73 = r\^l'^ j^^^ j The time for the reverse shock to cross 
the shell in the fixed frame is 

= (2) 



where a = 1/2 for RRS and a = 3/ Vl4 for NRS (Sari & Piran 1995). The shock expands to 
the radius i?A ~ cti\ — QiVTA when the reverse shock is just over. In addition to i?A there 
are two other characteristic radii involved here: the deceleration radius i?^ = ^?7~^, at which 
the swept-up mass is Mo/77, and the shell spreading radius Rs — ^rf. They have different 
dependences on A and 77 for the wind and ISM cases (Sari & Piran 1995). A parameter 



C — V introduced to compare these three radii. According to equation (1), we get 

Rr, = CRja = C'^., (3) 

with / = C'^rf. C 1 describes RRS: f "rf, and i?,, ^ _Ra <^ -Rs, while C ^ 1 describes 
NRS: / ^ rf, and Rg <^ R^ <^ i?^. For RRS, the previous assumption of A =const is 
always satisfied, and deceleration of the shell happens before the reverse shock crosses it. 
For NRS we must consider the spreading effect of the shell, and ^ — r/rf with 714 — Al/r^. 
The comoving density ratio f — l/r decreases with increasing r. The reverse shock becomes 
trans-relativistic at Rjsf — l/rf'{—R,j) when f — rf. K similar correction on and R/\^ must 
be done since A and / vary with r. We are most interested in RRS because of both a large 
energy output and a constant width. 

Before RRS crosses the shell, the distance that the shell travels when RRS crosses a 
length dx in the unshocked shell frame can be given by (Kobayashi 2000), 

dr — a'q^JJdx. (4) 

The shocked electron number is proportional to x/A. In the observer's frame we have r = 
2c7|i0, so the hydrodynamic variables are 



8^2 A /te\"^ SAuipC^ ft, 



-2 



where A^o — Mo/rUp is the total electron number of the shell and T — A/ (2c) is the observer's 
time at which the RRS crosses the shell, which is one half of the crossing time adopted by 
Kobayashi (2000). 

After RRS crosses the shell, the temperature of the shocked shell is very high and it will 
expands adiabatically and enters the spreading stage. Since it locates not far away from the 
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forward shocked wind, we take the BM self-similar adiabatic solution with impulsive energy 
injection for the evolution of the shocked shell (Kabayashi & Sari 2000). With ps oc ti^^^ , we 
get p3 oc oc r~^,73 oc r^^^-'^)/'^ oc r~^/^,n3 oc r^^*^"^^)/^ oc r~^/^ and r oc i®^. The 

hydrodynamic variables are therefore given by 

-3 = n3(T)(|)"'', e3 = e3(T)(|)"'', (7) 

// \ -3/8 

73 = 73(T) tI , iV, = iVo. (8) 



T 

The hydrodynamic evolution of a relativistic forward shock in a wind has already been 
discussed in some details (Dai & Lu 1998; Chevaher & Li 2000; Panaitescu & Kumar 2000). 
Before RRS crosses the shell, the Lorentz factor of forward shock ^sh — •\/272 is a constant 
and equals \/273. The forward shock radius is r = 2'y gj^ct^ = 47|ct0. Note that there is a 
difference of a factor of 2 between the reverse shock radius and the forward shock radius, 
which arises from the fact that the forward shock moves with 7^/1 at any time although the 
photon emitter travels with 72 (Chevalier & Li 2000). So the time when the reverse shock is 
finished becomes T' = T/2. For < T' 

8V2A /te\-' _8Amp^/t^\-^ 

^2 73/4 A. "SM \ rni I ' ^2 71 /OA 3/2 \ rri I 1 \^) 



^3/4^5/4 \rp, j > ;l/2A3/2 \T' 

with 72 = -\/2(Z/A)^/^/2. For > T", the hydrodynamic variables evolve as 

-2 = n2(T') (^)"'' , e2 = e2(T') (^)"'' , 72 = 72(T') (^)"'' . (10) 

The above analytical hydrodynamic evolution for both the reverse shock and forward shock 
is valid only in the relativistic phase. However, a late-time afterglow may be produced by 
a non-relativistic forward shock (Huang, Dai & Lu 1998, 1999, 2000). Thus, our results are 
applied only for early or very early afterglows in wind environments, including the reverse 
shock emission. 



3 LIGHT CURVE OF VERY EARLY AFTERGLOW 
3.1 Radiation from RRS 

The shock accelerates electrons to a power-law energy distribution with an index p and with 
a minimum Lorentz factor 7^ = Ce(73 — l)[(p ~ 2)/(p — l)]mp/me, where is a fraction 
of the internal energy density carried by the electrons. The standard radiation mechanism 
for GRB afterglows is synchrotron radiation by relativistic electrons in the magnetic field 
whose energy density is another fraction [^b) of the internal energy density. We also consider 
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synchrotron self-absorbtion at the radio band. Thus, the spectrum as discussed by Sari, Piran 
& Narayan (1998) in the ISM case can be used to discuss an afterglow in the wind case. 
The spectrum is characterized by three break frequencies: the typical frequency 1/^, the 
cooling frequency Uc, and the self-absorption frequency Va- A detailed derivation of Ua for 
fast-cooling and slow-cooling regimes is given in the appendix. Since B oc (^5/11)^/^73 and 
7c oc 1/ {-f^BH^), we have Pm oc 737™-^ and z/^ oc ^zllB. The peak flux F^^rnax oc Ne'ysB/D'^, 
where D is the distance from the GRB site to the observer. 
Before RRS crosses the shell (t® < T), 

= 8.65 X l^''E;,'''eedBUvlmA.i^'.^'" (I) Hz, (11) 
= 1.32 X 10-£;,f G?/2^-^aV^| Hz, (12) 

F-_, = 1.42E52?B(-2%oo^y' Ar3^D2t Jy, (13) 

where E = lO^^^Jss ergs, 77 = 3OO77300, A = lO^^Aia cm, = l^~^iB,-2, ie = O.GCe.o, 
^ = 3 X 10^^^* cm~^ and D — 10^^£)28 cm. The time for the reverse shock to cross the shell 
is T = 167Ai3 s. For i^c < i^a < ^^m, we have 

= 4.11 X l^^^Eli^^'Al'^^-,!'' (I) Hz, (14) 

which is consistent with eq. (18) of Dai & Lu (2001). For < i^c < we have 

= 6.47 X l^''E^r'%^'AT'^'.r''\t-. (I) " Hz. (15) 

The coohng Lorentz factor (7c oc t^/T) is determined by the dynamical timescale. At very 
early times, 7c may be less than unity, which is impossible in physics. In fact the radiation 
power of such electrons with 7c ^ 1 will be cyclotron power, instead of synchrotron power. 
According to Dai & Lu (2001), we neglect the radiation before t©,crit (same as to in Dai & 
Lu 2001) when 7c = 1, 

ie,crit = n.3E^2^'Ali'^B,-2A'J' s (16) 

which is the beginning for our calculation. 

When RRS crosses the shell {t^ > T), the shock-accelerated electrons will be in the 
slow coohng regime if Vc{T) > VmiT). However, there will be no radiation if Uc{T) < Vm{T), 
because all the previously shocked electrons have cooled and no newly-shocked electrons are 
in the shell when > T. The break frequencies and peak flux for i^dT) > i^m{T) are 

l^m = 1^7n{T) J ,l^cut^ ) ' ^"'"^^ ^ F^^rnax{T) J . (17) 
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We have also derived the self-absorption frequency: in the slow cooling phase, Ua oc {to^/T) 
for Ua<'^m< T^cuu and Va OC (tg5/T)-(i5p+26)/8(p+4) fo^, i/^ < z/„ < v^^f The time for the 
reverse shocked shell to become non-relativistic, i.e., when 73 = 2, is 

te,nr = m.^El^'A-^I^A-^^'^'T. (18) 
3.2 Radiation from the forward shock 

Before RRS crosses the shell, the forward shocked wind matter moves with the same Lorentz 
factor as the reverse shocked shell matter. Whilst most of the initial kinetic energy of the 
shell has been transferred to the forward shocked matter, as shown in eq. (3). The forward 
shock decelerates at a slower rate than the shocked shell does after the crossing time. We 
can consider the shocked shell as a tail of the forward shock based on the BM solution with 
its self-similar parameter x r^~^ oc r^, which describes inward motion of the shocked shell 
relative to the shock front as the shock radius increases. 

Before RRS crosses the shell (i© < T"), the radiation from the forward shock can be 
described by 

= 1.70 X l^^'EliX^eiU^-.!" (^) " Hz, (19) 
= 5.28 X 10-4f Cb^-^^^^A}/^^ Hz, (20) 

Fl^ma. = ^.QlEli^A^i^-^'^D-i Jy. (21) 

For t'c < < i^m, we have the self- absorption frequency 

/J- X -2/3 

= 2.92 X 10^^Ay3Ar3^/=^ ( I ) Hz, (22) 



and for t'a < t'c < we have 



= 9.78 X lQ^'E;^"A^^I'i^-l"et2 (^) " Hz. (23) 
The time for 7c > 1 is 

t'®,cru = 2.^lE^,'"i^'l,%,_,Al" s, (24) 

which is the beginning of our calculation. Since the peak flux of the forward shock emission 
is comparable to that of the RRS emission, a very early afterglow may be dominated by 
both of the two regions. An optical flash may even be dominated by the RRS emission if E 
is very large and A^ is small, which can be estimated by comparing eq. (13) with eq. (21). 
The break frequencies and peak flux for > T' are 
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The self-absorption frequency is: (1) in the fast coohng phase, Va oc (te/T')"^/^ for t'c < '^a < 
i^TO, and Va oc {t^/T')'^/^ for t'a < i^c ^ i^m! (2) in the slow cooUng phase, Va {t^/T')~^^^ 
for i^a < i^m < ^c- The case of v^ < Va < Vc seems impossible before ^0^^ = 2.763 x 

3.3 Light curves of optical and radio afterglows 

The above subsections give the spectral profiles and the decaying laws of the break frequen- 
cies. In this subsection, we will discuss the light curves at some observed frequencies such 
as Uopt — 4: X 10^^ Hz and Urad — 8.46 GHz. We calculate the evolution of the flux densities 
at these two frequencies with fixed parameters of ^^,0 — 1-0, Cb,-2 — 1-0 and A13 = 0.5 but 
with £'52 varying between 1.0 and 10^, and varying between 10~^ and 1, to investigate 
the effects of energy and wind intensity. The values of ,^e,o and ^^,-2 ai'c consistent with the 
afterglow of GRB 970508 (Wijers & Galama 1999; Granot, Piran & Sari 1999) and the opti- 
cal flash and afterglow of GRB 990123 (Wang, Dai & Lu 2000), and A13 = 0.5 is consistent 
with T' = A/4c ^ 42 sec, which is the peak time of the optical flash of GRB 990123. The 
light curves in six cases and at two frequencies are shown in Figures 1-6. 

The emission from RRS is nearly comparable or dominant over that from the forward 
shock at optical frequency. This is naturally estimated from the ratio of the reversely shocked 
electron number of the shell to the forward shocked electron number, 

^ = 3/V^ = 4.21£;,f .1-/2^-/^^3-1 , (26) 

before RRS crosses the shell, compared with eq.(58) of Ghevalier & Li (2000). In detail, the 
maximum optical flux in RRS reaches at — T, 

jprs lrp\ TTTS 

i/optV / v,max\ ^ J 

= 8.16 X 10-^£;rr;3-oU;^/'Ar3^/^eB:^2^2tJy, (27) 
when RRS is still in the fast cooling regime, as emphasized by Ghevalier & Li (2000). The 
above equation is not suitable for the case of E^2 = 10, A* = 10~^, where RRS enters the 
slow cooling regime before it crosses the shell, because we can use eqs.(ll) and (12) to get 
the transition time by = ^c, 

j-rs 

= 2.56 X 10'e;,'/\ooAI/'A^,'/'^,^o^s,-2 = 0.11 < 1, (28) 
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for = 10, A* = 10~^. The optical flux contributed by the forward shock reaches its 
maximum at — T', 

= 0.583EU'A^,'/%%tD^iJy, (29) 

except in the two cases of £'52 = 10, — 0.1 and £'52 = 10, A* = 10~^. For these two cases, 
the optical flux by the forward shock reaches the maximum before RRS has crossed the shell 
and decays to 

I'opt^ ' i>,raax\ / 

= l-583£,f (^)-/^(i^)n^/^_,D-Jy. (30) 

For more energetic flreball and weaker stellar wind, the emission from RRS becomes more 
stronger than that from the forward shock. The weakest RRS case is for E = 10^^ ergs, 
A^. = 1.0, where the forward shock dominates RRS at Uopt by about one order of magnitude 
(Fig. 2.). The peak flux F^^^^^p^ak increases with both E and A^ and varies from ~ 1 Jy to 
~ 10^ Jy. An optical flash occurs in all cases, although its properties such as the brightening 
and decaying timescales and the light curve index may be different from each other. The 
time scale of emission with flux above 1 Jy ranges from 10 to 10^ seconds. None of the six 
cases discussed here is consistent with the rapid brightening {F^^^^ oc t^'^) and fast decaying 
(oc t^^'^) of the optical flash of GRB 990123. This further conflrms the result of Chevalier 
& Li (2000) that the optical flash of GRB 990123 originates from a homogeneous ISM. 

The radio afterglow is mainly attributed to radiation from the forward shock and we can 
neglect the contribution of the emission from RRS. The light curve indices of the forward 
shock are 2, —2/3, 0, — (3p — l)/4 from low- to high-frequencies in Figs. 1, 4 and 6, as 
shown in Chevaher & Li (2000). We also get an index P — 3, which does not equal to 7/4 
of Chevalier & Li (2000) at very early times. Indices of 2, 1, and in Fig. 2. and Fig. 5. are 
consistent with the results of Chevalier & Li (2000). The peak flux at radio band varies by 
two orders of magnitude, from 10^ to 10^ /iJy. There is a long lasting platform in the radio 
flux before the forward shock enters the non-relativistic phase ( also see Fig. 3 of Panaitescu 
& Kumar 2000). It might be used to distinguish the wind environment from the uniform 
ISM environment of GRBs and give further hints on the central engine. 
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4 SUMMARY 

We have analyzed very early afterglows in wind environments by considering the radiation 
from both a relativistic reverse shock and a relativistic forward shock. The resulting optical 
flash in the wind case is mainly attributed to RRS because of reduction of the flux from the 
forward shock by synchrotron self-absorption. Fo more energetic fireball and weaker stellar 
wind, the optical emission from RRS becomes more stronger than that from the forward 
shock. The resulting radio flare is largely produced by the forward shock emission. 

Theoretically, the peak flux of the optical flash ranges from 1 Jy to 10^ Jy, and the 
peak time since the GRB trigger is tens of seconds, which is determined by the shell width. 
In our cases we choose the width according to the typical time scales of long gamma-ray 
bursts, i.e., ~ a few tens of seconds. Such a large optical fiux at very early times is due to 
the dense wind at small radius. If optical observations could be done in a few seconds after 
the GRB trigger, then one can discriminate the wind environment from the uniform ISM 
environment. It is a challenge for current observational instruments. While strong optical 
flashes have been predicted theoretically in this research as well as by many other authors, 
in realistic observations optical flash has been observed only from GRB 990123 till now. The 
missing of a significant number of strong optical fiashes has been interpreted as due to dust 
obscuration (Soderberg & Ramirez- Ruiz 2002). We propose that further contribution to the 
reduction of the optical flash may come from the dense shell around the stellar wind with the 
column density about 10^° — 10^^ cm~^ (Ramirez- Ruiz et al. 2001). Another characteristic of 
wind environment which might be easier for observations is a long lasting platform in radio 
flux, which arises from the interaction of the forward shock with the wind matter. 

However, dust extinction could not darken a 4th magnitude (~ 100 Jy) optical flash sig- 
nificantly to be immune to the optical survey. The non-detections of optical fiashes brighter 
than ~ 9th magnitude would eventually constrain the isotropic energies of GRBs to be no 
more than a few 10^^ ergs and result in relatively weak circum-burst wind intensities with 
A* < 1.0, according to eqs. (13), (21) and (27). We expect future observations on early 
afterglows to diagnose the environments of GRBs efficiently. 
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APPENDIX A: SYNCHROTRON SELF- ABSORPTION BY ELECTRONS 
WITH A BROKEN POWER-LAW DISTRIBUTION 

Al Synchrotron self-absorption coefficient by electrons with one single 
power-law distribution 

All the quantities below are in the comoving reference frame. The power law distribution of 
the shock-accelerated electrons can be written as N{'ye)d^e — A^TeTj^^Te, where N-y^ is the 
normalized coefficient and 71 < 7e < 72- The power radiated from an electron with 7e is 
(Rybicki & Lightman 1979) 

P(^,7e)^ ^ / ^5/3(^)0?^, (Al) 

where ul — QeB /2Tr'meC is the Lamor frequency and Uc = (37gZ/2, sin^)/2 is the typical 
frequency emitted by this electron. We rewrite the self-absorption coefficient of eq. (6.52) of 
Rybicki & Lightman (1979), and get 



STTTTleiy^ J -11 d'ye \ 7e / 

Before integrating the above equation, we define ui — {SjIul sin 9)/2 and 1/2 — {2>^Ivl sin 0)/2. 
For i/ <S z/i, we have 

For simphcity, we take an isotropic distribution of 9 with < (sin^)^/^ >= '^Z^^} and obtain 

5i (,g j 

K = 114.65 X ^7V^^s2/37-(^'+2/3)^-5/3_ (A4) 
For i/i ^ 1/ <^ 1^2, we obtain 

= (">.c^)'-'iV,.i3<'+^)/^(sinO)<'«)/V-<'«)/2, (A5) 

where g{j)) = ^V{^^)V{ ^^^^^ ){j) + y). This case is the same as the well known eq. (6.53) 
of Rybicki & Lightman (1979). For V2 ^ ^, we get 
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and if we take the isotropic distribution < sin^^ >= 1/4, 



9 . / re ,z ^ 

which is always neglected in astrophysics but it may be applied for GRBs in which the self- 
absorption i/a > ^2 is possible and it may explain the spectrum of GRBs combining inverse 
Compton scattering. 

For unification and convenience for comparison in the next section we rewrite eqs. (A3), 
(A5) and (A7), 

K = CiIat^^Ti"^^'^ i-V" ,^^^1 (A8) 



ku = C2|iv,^7r^'^'^ (-) .yi<^y « (A9) 



K = C3|Ar^,72"^"^'^ ( - ) e-^l^\ « V, (AlO) 



3 



where c\ = — '1^""^ ^ , C2 = ^22g[p) and C3 = ^^"^^ (j9 + 2) which are of the same order. 
For illustration, ci = 14.78, C2 = 17.80 and C3 = 13.64 for p = 2.5; Ci = 15.60, C2 = 15.64 and 
C3 = 12.12 for p = 2.0. So the above expression of k^, still is valid even when u approaches ui 
or 1/2. For simulations, the discontinuity of can be removed if we constrain ci — C2 — cs/e. 

A2 Self-absorption frequency and the spectrums of fast and slow cooling cases 

Now we consider the absorption by electrons with a broken power-law distribution. The 
break Lorentz factors 71 < 72 < 73 with power law indices pi of (71,72) (hereafter we 
refer to this region as region A) and p2 of (72,73) (region B). For slow cooling pi — —p 
and P2 — —{p + 1) with 71 = 7^, 72 = 7c, and 73 = 7m- For fast cooling pi — —2 and 
P2 — — (p + 1) with 7i = 7c, 72 = 7m, and 73 = 7m- Both regions A and B contribute to k^- 
According to eqs. (A8), (A9) and (AlO), we find that ki, is mainly determined by a nearer 
region with respect to its frequency u, i.e., k,y = ki,{A) for i' < I'l and Ui < u < 1^2, while 
ku = ky{B) for z/2 < z/ < z/3 and < u. 

The self- absorption frequency Ua is determined by /c^/^L 1, where L is the length of the 
radiation region in the comoving frame. The spectrum below Ua will be corrected since the 
optical depth is larger than unity and I^, oc j^/k^. From the above analysis, we obtain the 
complete spectrum in all cases. 

In the fast cooling phase, (1) for Ua < i^c < ym the spectral indices of four segments 
are (2, |,— |,— |); (2) for Vc < < ym the spectral indices are (2, |,— |,— |); (3) for 
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J^c < '^m < ^a(< J^m) the spectral indices are (2, |,|,— |). Case (3) is previously rarely 
discussed in astrophysics and maybe contributes to GRBs but not to their afterglows. There 
is still case (4) um < Va- We do not discuss it since the flux oIv>vm is negligible. 

In the slow cooling phase, (1) for < Vm < i'c{< ^m) the spectral indices are (2, |, — — |); 
(2) for i^m < i^a < T^c{< t^m) the spectral indices are (2, |, — — f ); (3) for v^n < Vc < Va{< 
Um) the spectral indices are (2, |,|,— |). For case (1), from eq. (A4) we can flnd Ua in 
the observer's frame is consistent with eq. (22) of Wijers & Galama (1999), or we use cq. 
(A8) to get the same z/^ as eq. (52) {u < i/p) of Panaitescu & Kumar (2000), if we take 
N^^ ~ (p — l)ne7^~^ s-^id is the electron density in the comoving frame. 
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Figure 1. Light curves of very early afterglows in a stellar wind. The upper panel corresponds to optical flashes at i/opt = 4x 10 
Hz, and the bottom panel to radio flares at 8.46 GHz. The dashed, dotted and solid lines represent the radiation from the 
relativistic reverse shock and forward shock and their total emission. E = lO^^ergs, A = 3 x 10^* cm~i , rj = 300, A = 5 x lO^^cm, 
= 0.6, = 10~^, D = 10^* cm and p = 2.5 axe assumed, tea ^ 856sec is the time when Vc = Vopt- 
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Figure 3. Light curves of very early afterglows in a stellar wind. Same as Fig.l except for E = lO^^ergs, A = 3 x 10^^ cm ^ . 
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